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Abstract 

A finite set A in a complex sphere is called an s-code if the number of inner products 
of two distinct vectors in X is equal to s. One of major problems on s-codes is to classify 
the largest possible 3-codes for a given dimension. A complex spherical 3-code naturally 
has the structure of an oriented graph. In this paper, we study the minimal dimension 
which contains a complex 3-code having the structure of a given oriented graph. It is 
possible for certain oriented graphs to determine the minimal dimension. The algorithm 
to classify the largest 3-codes can be obtained by using this minimal dimension. The 
largest 3-codes for dimensions I, 2, 3 are classified by a computer calculation based on 
the algorithm. Roy and Suda (2014) gave some upper bounds for the cardinality of a 
3-code, and the set attaining the bound is said to be tight. A tight code has the structure 
of a commutative non-symmetric association scheme. We classify tight 3-codes for any 
dimension. Indeed there exists no tight 3-code except for dimensions 1,2. 

Key words: complex spherical s-code, s-distance set, tight design, extremal set theory, graph 
representation, association scheme. 

1 Introduction 

Let A be a finite set in the d-dimensional complex unit sphere U(d). The angle set D(X) is 
dehned to be 

D{X) = {x*y \ x,y G X,x ^ y}, 

where x* is the transpose conjugate of a column vector x. A finite set X is called a complex 
spherical s-code if \D{X)\ = s and D{X) contains an imaginary number. The value s is called 
the degree of X. For A, X' C U(d), we say that A is isomorphic to X' if there exists a unitary 
transformation from A to A'. An s-code A C U(d) is said to be largest if A has the largest 
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possible cardinality in all s-codes in Q(d). One of main problems on s-codes is to classify the 
largest s-codes for given s and d. 

An s-code in a real sphere S'^~^ is known as an s-distance set, that has s Euclidean 
distances. For an s-distance set X in Delsarte, Goethals, and Seidel [B] gave an upper 

bound 

An s-distance set X is said to be tight if X attains this bound. A tight s-distance set has the 
structure of a Q-polynomial association scheme, and becomes a tight spherical 2s-design [B]. 
Tight s-distance sets have been classified except for s = 2 [iiEiiaE]. The largest 1-distance 
set in is the regular simplex. The largest s-distance set in 5^ is the regular (2s -|- 1)- 
gon. The largest 2-distance set in 5"^“^ has been determined for d < 93,d jtz 46^ 78 BUM- 
The largest 3-distance set in has been determined for d = 3, 8, 22 [THl HD]. The largest 
spherical s-distance set is not known for other (s,d). The classification of largest spherical 
s-distance sets is still open except for (s,d) = (l,d), (5,2), (2,d < 7), (2,23), (3,3). 

For a 2-code X in n(d), we have an upper bound 

2d -|- 1 if d is odd, 

2d if d is even. 



A 2-code X is said to be tight if X attains this bound. For odd d (resp. even d), the existence 
of a tight 2-code in Q(d) is equivalent to that of a doubly regular tournament (resp. skew 
Hadamard matrix) of order d [T5|. For a 3-code X in n(d), we have an upper bound 

, , fd if d = 1, 

1^1 < S 9 

l^d^ -|- 2d if d > 2. 

A 3-code X is said to be tight if X attains this bound. Roy and Suda m proved that a tight 
3-code has the structure of a commutative non-symmetric association scheme. In this paper, 
we classify the tight 3-codes in Q(d) for any d. Indeed there exists no tight 3-code except for 
d = 1,2. Moreover we classify the largest 3-codes in Q{d) for d = 1,2,3. 

We use complex representations of oriented graphs in order to classify the largest 3-codes 
in 11(d). An oriented graph is a directed graph which has no symmetric pair of directed edges. 
An oriented graph G = (F, E) is representable in 11(d) if there exist a mapping (p from V to 
11(d), an imaginary number a with Im(Q!) > 0, and a real number j3 such that for all u,v gV, 


p{u)*(p{v) 


a if {u, v) G E, 
< a if (u, u) G E, 
/3 otherwise. 


The image of this map p is called a eomplex spherieal representation of G. For two oriented 
graphs G and G', if G is not isomorphic to G', then a representation of G is not isomorphic to 
that of G'. The Gram matrix id of a complex spherical representation of G can be expressed 
by 

H = M + cV^{A- A^), 

where c is a real number, A is the adjacency matrix of G, and M is a real matrix. This matrix 
M is the Gram matrix of a real spherical representation of a simple graph G' with adjacency 
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matrix A + A'^. To make the theory of complex representations of oriented graphs, we will 
recall known results of Euclidean representations of a simple graph in Section [21 and prove 
several results on Hermitian matrices in Section [3l We determine the smallest dimension of 
the complex representations of an oriented graph satisfying a certain property in Section 31 
An algorithm to classify the largest 3-codes is given in Section [5l and the largest 3-codes in 
Q{d) is classihed for d = 1,2,3 by computer calculation. In Section [6l we classify the tight 
3-codes in Q(d) for any d by the association schemes attached to the codes. 

2 Euclidean representations of a simple graph 

In this section, we introduce several results about the embedding into of a simple graph. 
They are used later to give a representation of an oriented graph. Let G = (V, E) be a simple 
graph which is neither a complete graph nor isolated vertices, where \ V\ = n. Let A be the 
adjacency matrix of G, and A that of the complement. The matrix Mc{G) = Me is defined 
to be 

Me{G) =cA + A, 

where 0 < c < 1. A graph G is said to be representable in if Me is the distance matrix of 
some hnite set X in for some c. This embedding X is called a Euclidean representation of 
G. Let Rep(G) be the smallest d such that G is representable in M'^. 

Theorem 1 ([7]). For a simple graph G, there exists 0 < ^ < 1 such that the following hold: 

(1) Mg is the distance matrix in Rep(G) dimension. 

(2) For ^ < c < 1, Me is the distance matrix in n — 1 dimmension, and not in n — 2 
dimmension. 

(3) For 0 < c < Me is not a distance matrix in any dimension. 

The embedding Mg in Theorem [H is called the minimal representation of G. Roy [18] 
determined Rep(G) as follows. Let j be the all-ones column vector. 

Theorem 2 f |18[ Lemmas 4,5,6, Theorem 7]). Let Ti he the i-th smallest distinct eigenvalue 
of A, rui the multiplicity of Ti, and £i the eigenspace of Ti. Let Pi be the orthogonal projection 
matrix onto Si. Let fdi be the main angle of Ti, namely, Pi = a /{Pi ■ {Pi ■ j)/n. Then the 
following hold: 

(1) If Pi = 0, then C = ('Ti + ^)/'Ti and Rep(G) = n — mi — 1. 

(2) If Pi 0 and mi > 1, then .^ = (ti -|- 1)/ti and Rep(G) = n — mi. 

(3) If P 2 = 0, mi = 1, T 2 < -1, and PI/{t 2 - n) = X]i >3 Pf/ip - T 2 ), then f, 
and Rep(G) = n — m 2 — 2. 

(4) If P 2 = 0, mi = 1, T 2 < -I, and PI/{t 2 - n) > 'P2i>3 C 

and Rep(G) = n — m 2 — 1. 

(5) Otherwise, we have ^ < (ri -|- l)/ri, {t 2 + 1)/7‘2 and Rep(G') = n — 2. 

A graph G is said to be of Type (i) if G satishes the condition (i) in Theorem |2] for 
i = 1,... ,5. A representation of G is said to be spherical if it can be on a sphere. 


= (t 2 -k 1 )/ t 2 
= {t2 + 1)/t2 
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Theorem 3 f|13jh Let G be a simple graph. Then the following hold: 

(1) If G is of Type (1), (2), or (4), then the minimal representation of G is spherical. 

(2) If G is of Type (3) or (5), then the minimal representation of G is not spherical. 

(3) Representations in Theorem [T] (2) are spherical. 

A non-negative symmetric nxn matrix with zeros on the diagonal is called a dissimilarity 
matrix. The embedding dimension of a distance matrix M is the smallest dimension d, such 
that there exists a set in of the distance matrix M. Let J denote the matrix whose entries 
are all one. 

Lemma 1 ([E])- Let M be a dissimilarity matrix. Then the following are equivalent. 

(1) There uniquely exists a > 0 such that —M + aJ is a positive semidefinite matrix of rank 
d, —M + a'J is a positive semidefinite matrix of rank d+1 for any a' > a, and —M + cJ 
is not positive semidefinite for any c < a. 

(2) M is the distance matrix of a set in where d is the embedding dimension of M. 

3 Results on Hermitian matrices 

In this section, we prove several results on Hermitian matrices. Let iL be an n x n Hermitian 
matrix with s distinct eigenvalues ri < T 2 < • • • < r^. Let Si be the eigenspace corresponding 
to Ti, and Pi the orthogonal projection matrix onto Si. The main angle fli of Tj is defined to 
be /3j = {Pi ■ j)*{Pi ■ j)/n. Note that /3i = 0 if and only if Si C j^. An eigenvalue r* is said 
to be main if /3j ^ 0. Let be the multiplicity of Ti. Let P denote the orthogonal projection 
matrix onto j'^, namely P = I — (1/n) J, where I is the identity matrix. 

Theorem 4 1[15]L Let H be a Hermitian matrix of size n, and M = H + aJ, where a is a 
real number. Let ti < T 2 <■■■< Tr be the distinct main eigenvalues of H, and Pi the main 
angle of Ti. Let p-i < fi 2 < ■ ■ ■ < Rs be the distinct main eigenvalues of M. Then r = s holds, 
and 

r r r ^q2 

Ylifii-x) = '[l{Ti-x){l + aJ2:p-^)- (1) 

i=l i=l j=l ^ 

Moreover, if a > 0, then ti < < T 2 < ■ ■ ■ < Tr < p.r, o,nd if a < 0, then p.i < ti < fj ,2 < 

■■■< Hr < Tr- 

Lemma 2. Let H be a Hermitian matrix. Let ti < T 2 < ■ ■ ■ < Tr be the main eigenvalues of 
H. Suppose H is not positive semidefinite. Then the following are equivalent. 

(1) There exists a > 0 such that H + aJ is positive semidefinite. 

(2) T 2 > 0, Y7i=i I Ti < 0, and PHP is positive semidefinite. 

Moreover, if {!) holds, then a > — l/(Ei=i holds. 
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Proof. Let r be an eigenvalue of H with main angle zero, and v an normalized eigenvector 
corresponding to r. 

(1) ^ (2): Since r is not a main eigenvalue, v is orthogonal to the all-ones vector. By the 

assumption of (1) that H+aJ is positive semidefinite, we have r = v*Hv = v*P{H+aJ)Pv > 
0. Since H is not positive semidefinite, we have ri < 0. Let < /i 2 < • • • < /^r be 
the main eigenvalues of if + aJ. By Theorem HJ we have ti < /ii < r 2 . Since H + aJ 
is positive semidefinite, we have 0 < /ri. Thus T 2 > 0, and by ([I]), < 0 

a > /p)- (In particular, /ii = 0 if and only if a = n/3j^/rj) > 0.) 

Since H + aJ is positive semidefinite, so is P{H + aJ)P = PHP. 

(2) ^ (1): Since v is orthogonal to the all-ones vector and PHP is positive semidefinite, 
we have 


r = v*Hv = v*PHPv > 0. (2) 

Since H is not positive semidefinite, we have ri < 0 and mi = 1. By m and T 2 > 0, a matrix 
iL -I- oJ is positive semidefinite for a > — l/(X]i=i >0. □ 

We can verify the following remarks by the proof of Lemma [2j 

Remark 1. If Lemma[2] (1) holds, then 

(1) Rank(R -|- aJ) = Rank(R) — 1 for o = —l/(X]i=i /p)^ 

(2) Rank(R -I- aJ) = Rank(R) for a > — l/(X]i=i /p)- 

Remark 2. If Lemma [2] (2) holds, then the null space of H is contained in j-^. 

Remarks. If Lemma[2](2) holds, then Rank(R-l-aJ) = Rank(PRP) fora = — l/(X]i=i 

A Hermitian matrix H can be expressed by R = M -|- y/—lA, where all entries of M and 
A are real. Let £q be the null space of M, and Tq that of \/—\A. We use this notation in the 
following theorems. 

Theorem 5. Suppose H is a positive semidefinite Hermitian matrix. Then £q C £o holds. 

Proof. Since M is a real symmetric matrix, we can take a basis of £q consisting of real vectors. 
For real v G £q, 

v*Hv = v*Mv + \/—lv*Av = 0, 

because A is skew-symmetric. Since R is a positive semidefinite, v*Hv = 0 if and only if 
Hv = 0. Thus we have 

0 = Hv = Mv -|- y/—lAv = i/—lAv. 

Therefore £'q C Tq holds. □ 

Theorem 6. Suppose H is a positive semidefinite Hermitian matrix. Then 2Rank(R) > 
Rank(M). 

Proof. By Theorem^ we have £q C Sq. If there exists non-zero v G {£Q)-^r\£o, then v*Hv > 0. 
Let T+ be the direct sum of eigenspaces corresponding to the positive eigenvalues of y/—lA, 
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£- that corresponding to the negative eigenvalues. It is easily proved that dim£l_|_ = dim£l_. 
For u G £1+ © {{£q)^ n £q), we have v*Hv > 0 because M is positive semidefinite. Therefore, 

Rank(FI) > dim(T+ © ((Tq)'*' Fl To)) 

= dim(T+) + dim((To)-^ n To) 

> — Rank(T) + (Rank(M) — Rank(^)) 

= Rank(M) — — Rank(^) 

> Rank(M) — - Rank(M) 

= ^Rank(M). □ 

Theorem 7. Let H he a Hermitian matrix. Suppose M is positive semidefinite, and Tg C To 
holds. Then there uniquely exists tj > 0 such that the following hold: 

(1) M + rjyJ—lA is positive semidefinite, and its rank is smaller than Rank(M). 

(2) For 0 < c < T], M + c\/—lA is positive semidefinite, and its rank is equal to Rank(M). 

(3) For 7] < c, M + c\/—lA is not positive semidefinite. 

Proof. We define a function 

<k(c) := min v*{M + c^/—lA)v. 

,v*v=l 

Note that <k(c) > 0 if and only if M + c^/—lA is positive semidefinite, and Rank(M + 
c^J—lA) < Rank(M). In particular, <k(c) = 0 if and only if Rank(M + c^J—lA) < Rank(M). 
Since <k(c) is the minimum value of the collection of linear functions in c, <k(c) is convex. 
Since M is positive semidefinite, we have ‘k(O) > 0. There exists v G (Tg)"*" such that 
V* —lA)v < 0. Therefore, we have limc^.oo ‘h(c) < 0. By the intermediate value theorem, 
this theorem follows. □ 


4 Representations of an oriented graph 


A complex spherical 3-code X has an angle set D{X) = {a,a,/?}, where a is an imaginary 
number with Im(a) > 0, and /I G M. A 3-code X gives an oriented graph G = {X,E), 
where E = {{x,y) G A x A | x*y = a}. Let G' be the simple graph {X,E'), where 
E' = {{x,?/} I (x,y) G E}. The adjacency matrix A of a graph is the matrix indexed by the 
vertices, with entries 


A 


xy — 


I if {x,y) in the edge set, 
0 otherwise. 


Let A, B be the adjacency matrix of G, G', respectively. Let B be the adjacency matrix of 
the complement of G'. The Gram matrix H oi a complex spherical representation of G can 
be expressed by 


H = M + cy/^{A- A'^), 
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where c is a real number, and M is a real matrix. Let denote a map from Q{d) to 5^“^ ^ 
defined by 

+ 2/1 v^, ...,Xd + ydV^) = Xd, yd)- 

Note that (j){x)'^(j){y) = Re(x*y) for x,y & Q.{d). The matrix M is the Gram matrix of 
ip{X) = {ip{x) I X E X}, and a spherical representation of G'. By Lemma [U we can express 
for some a > 0, 6 > 0, 

M = -{hB + B) + aJ 

after rescaling the long distance to 1. Since M is spherical, M is the minimal representation 
of G' of Type (1), (2) or (4), or M is not a minimal representation by Theorem [3l 

By Theorem [5l the null space of M must be contained in that of y/—l{A — A^). This 
condition is very strong, and the minimal representation of G' hardly satisfies it. We give 
simple examples: 


( ° 

1 

0 

0 \ 


( ° 

0 

0 


0 

0 

1 

0 

, G2 : A2 — 

1 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

V 1 

0 

0 

0 y 


V 1 

0 

0 

0 / 


Then both G'l and G 2 are the cycle (74. Indeed is of Type (1), and its minimal represen¬ 
tation is the vertex set of the square in R^. The Gram matrix of the square can be expressed 
by 

( \ 0 0 \ 


Ml — + B) + - J 


0^0-. 
0 b o' 


V 0 -i 0 i / 

The null space of Mi is Span{(l, 0,1,0), (0,1, 0,1)}. This coincides with the null space of 
Actually we can give a minimal-dimensional representation in 14(1) of Gi as 


Hi — -{-^B + B) + - J -I- — AJ) 


( 1 

/ 2 2 

1 

2 2 

2 2 

\ _1 

^2 2 


1 

7=T 

\ 

2 

2 

V-I 

1 


2 

2 


1 



2 

2 


V-1 

9 

1 

2 

/ 


On the other hand, the eigenvalues of \/~l (^2 — A^) are and hence 

the null space is empty. In this case, Rank(M 2 ) must be 4, and we use a non-minimal 
representation of G': 


M2 


-{B + B) + J = 


/ 1 0 0 0 \ 
0 10 0 
0 0 10 
\ 0 0 0 1 / 


Then we can give a minimal-dimensional representation in 14(2) of A 2 as 


H 2 


( 


1 


-{B + B) + J + 



-1{A2-aI) = 



\ 





0 
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In this section, we give the minimal-dimensional representation of an oriented graph G 
obtained from the minimal representation of G'. The dimension of a representation of G with 
a non-minimal representation of G' is not so small by Theorem [6l Therefore it will be enough 
to consider representations obtained from the minimal representation of G' to determine the 
largest 3-codes. 

Throughout this section, we suppose G' has non-zero B and B, and is of Type (1), (2), 
or (4). For ^ > 0 given in Theorem [H and a,c > 0, the matrix H{a,c) is defined to be 

H{a,c) =-{CB+ B)+ aJ+ c^/^{A-A^). (3) 

We want to determine a > 0 and c > 0 such that H{a, c) is positive semidefinite, and the rank 
of H (a, c) is minimal. Note that the matrix ^B + B is the distance matrix of the spherical 
minimal representation of G'. Let Sq be the eigenspace corresponding to the zero eigenvalue 
of ^/—1{A — A"^), and £q be that of —{^B + B). 

Remark 4. If G' is of Type (1), (2), or (4), then £q C j-^ holds by LemmaHJand Remark[2j 

If H{a, c) is positive semidefinite, then H{0, c) satisfies Lemma[2](2), and PH{0, c)P is pos¬ 
itive semidefinite. By Remarks [T] and El if P[{0,c) satisfies Lemma El (2), then there uniquely 
exists a > 0 such that Rank(R(a, c)) is minimal, and Rank(R(a, c)) = Rank(PR(0, c)P). 
Therefore we would like to choose c such that PH (0, c)P is positive semidefinite, and Rank(PR(0, c)P) 
is smaller. The following lemma shows such possible c and the evaluation of Rank(PR(0, c)P). 

Lemma 3. Notation is defined as above. Suppose £q C £o. Then there uniquely exists rj > 0 
such that 

(1) PH{0,r])P is positive semidefinite, and 

Rank{PH{0,rj)P) < Rank(PR(0,0)P), 

(2) PH{0,c)P is positive semidefinite, and 

Rank(PR(0, c)P) = Rank(PR(0, 0)P) 

for any 0 < c < r], 

(3) PH{0, c)P is not positive semidefinite for any rj < c. 

Proof. We have 

PH{0, c)P = -P{(B + B)P + c\^P{A - A^)P. 

It is easily shown that the null space of —P{(,B+B)P is contained in that of y/—lP{A 
By Theorem [71 this lemma follows. 

Next we have to check whether H{0,c) satisfies Lemma [2l (2) for 0 < c < r/ in Lemma [3l 
If H{0,c) satisfies Lemma [2l (2), we can construct a representation of G by choosing suitable 
a. Let 

U = {(a, c) I H{a,c) is positive semidefinite, a > 0, c > 0}, 

and 

Rep(G) = min{Rank(R(a, c)) [ (a,c) G U}. 


-A^)P. 

□ 
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Theorem 8. Notation is defined as above. Suppose £q C Sq. Let rj he defined in Lemma^ 
Then the following hold. 

(1) If H{0,t]) satisfies Lemma [2] (2), then 

Rep(G) = Rank(iL(0, r/)) — 1 < Rep(G'). 

(2) If H[0,7]) does not satisfy Lemma [2] (2), then 

Rep(G) = Rank(iL(0,0)) - 1 = Rep(G'). 

Proof. (1): For H{0,c) satisfying Lemma [5] (2), there exists Oc > 0 such that H{ac,c) is a 
positive semidefinite matrix whose rank is Rank(Lf(0, c)) — 1. If iL(0, c) satisfies Lemma [2] 
(2), then it is satisfied that 0 < c < r/ by Lemma El Suppose there exists 0 < c < rj such that 
H(0,c) satisfies Lemma [2] (2). By Lemma [3] and Remark El we have 

Rank(Lf(ac, c)) = Rank(PLf(0, c)P) 

= Rank(Pi/(0, 0)P) = Rep(G') 

> Rank(PLf(0, ry)P) 

= Kank{H{orj, r])) 

= Rank(Lf(0, r/)) — 1. 


Therefore (1) holds. 

(2): Since G' is of Type (1), (2), or (4), Lf(0,0) satisfies Lemma[21(2). By Remarkdl there 
exists a' > 0 such that II(a', 0) is a positive semidefinite matrix whose rank is Rank(Lf(0,0)). 
Since Lg C j-^ by Remark El the null space of Lf(a', 0) is also Lg. By Theorem [71 there exists 
T]' > 0 such that H{a',r]') is positive semidefinite. This implies that there exists 0 < t]' < rj 
such that Lf(0,r/') satisfies Lemma[2l(2). For such r/', we have 

Rank(R'(a^/,r?')) = Rank(P77(0, r/')R) = Rank(P77(0, 0)P) = Rep(G') = Rank(R'(0,0)) - 1 
by Remarks mEl and LemmaEl(2). Therefore (2) holds. □ 

5 Algorithm to give the largest 3-codes 

In this section, we classify the largest complex 3-codes in 11(d) for d = 1, 2,3 by the following 
algorithm. First we classify the simple graphs G' = (X,E') of Type (1), (2), or (4) with 
Rep(G') < 2d by the algorithm in [5]. Indeed the classihcation is possible when Rep(G') < 8 
by a current computer. For such a graph G', we give the distance matrix f^B+B of the minimal 
representation, where B is the adjacency matrix of G', and B is that of the complement G'. 
The null space Lg of — + R) is determined as follows: 


II 

if G' is of Type (1) 

H 

C 

II 

u? 

if G' is of Type (2) 

II 

if G' is of Type (4) 


by Theorem El We give oriented graphs G such that A + coincides with R or R, where 
A is the adjacency matrix of G. By Theorem El we would like to obtain all oriented graphs 
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G such that £q C £o, where £o is the null space of — A^). In particular, any row in 

A is perpendicular to any vector in £q. For G satisfying £q C £o, we can determine Rep(G) 
by Theorem [HI The condition £q C £o is very strong, and it reduces the huge number of 
possibilities of A. 

The Gram matrix of a Euclidean representation of + B has an expression —{^B + 
B) + aJ for a > —l/(X][=i by Lemma[2]and Remark[TJ Let £q be the null space of 

— {B + ^B) + aJ, where a = — l/(X]/'=i Note that £q C £q. For a simple graph G' 

with Rep(G') = 2d, we give oriented graphs G such that £q C £q. By Theorems [5] and (H for 
G' with Rep(G') = 2d, the condition Tq C Tq is enough to give a complex spherical 3-code in 
£l{d). Note that the condition £q C £o is much stronger than £q C Tq. 

We define the following sets of oriented graphs up to isomorphism. 

Bi{d) := {G \ £'o C £o,Rep{G') = d}, 

B 2 {d) := {G I £'' C To,Rep(G') = d}, 

2d-l 

B{d) :={Ge IJ U B 2 { 2 d) \ Rep(G) = d}, 

i=d 

and m := max{|R| | {V,E) G B{d)}. We used the function in Magma to check the isomor¬ 
phism of oriented graphs. In our algorithm, we deal with a representation of G obtained 
form the minimal representation of G' . Since we ignore representations of G obtained from 
non-minimal representations of G', if m < 2d + 1, then the classification of the largest 3- 
codes in £l{d) cannot be given with our algorithm by Theorem [6l If m > 2d + 1, then 
{{V,E) G E{d) I |F| = m} is the classification of largest 3-codes in Q{d). 

By the above algorithm we can classify the largest complex 3-codes in £l(d) for d = 1,2, 3. 
Table 1 shows the number of largest 3-codes. 


d 

1 

2 

3 

|A| 

4 

8 

9 

# 

1 

1 

35 

r 

"able 1 



For d > 4, a usual computer cannot give the classificaition. For d = 1, 2, the largest complex 
3-codes are tight, which will be given in Section [6j For d = 3, one of the largest 3-codes is the 
union of three equilateral triangles in C^, which are orthogonal to each other. For the other 
largest 3-codes X, (j){X U u is the unique largest 2-distance set in R® 

mm, which is the minimal representation of the Schlafli graph with 27 vertices. 

6 Largest complex spherical 3-codes 

In this section, we give bounds on complex spherical 3-codes and characterize the tight 3- 
codes. In m another type of codes, 5-codes, are defined, see m for more details. Here 
5 is a lower set of N^, meaning that 5 is a finite set, and if (f,j) G is in 5 then so is 
{k,l) for any Q < k < i and 0 < I < j. A finite set X in 0(d) is called an S-code if there 
exists a polynomial E{x) = Yl,{k i)&s with real coefficients such that E{a) = 0 for any 

a G D(X). Let Hom(/c, 1) denote the set of functions representable by polynomials which are 
homogeneous of degree k in variables {zi ,..., z^} and homogeneous of degree I in variables 
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{zi,... ,^} for z = (zi,..., Zd) G 0(d). Let Harm(fc, 1) denote the subspace of Hom(A:, 1) that 
is the kernel of the Laplace operator A = Yli=i d'^ldzidzi. An upper bound on the size of an 
5-code is given as follows. 

Theorem 9 ([T71 Theorem 4.2 (ii)]). For d > 2, let X be an S-code in 0(d). Then \X\ < 
E„,„, 5 dim(Harm(t,i)) = f'fe') {"Ih') - (TT") fth") 

An 5-code is said to be tight if equality holds. For a finite lower set T, X is called a 
complex spherical T-design if, for every polynomial / G Hom(A:, 1) such that {k, 1) is in T, 


1 

A 


E /(") = / 


f{z)dz, 


(4) 


where dz is the Haar measure on 0(d) normalized by = 1. As shown in the following 

theorem, tight 5-codes are complex spherical 5 * 5-designs, where U *V := {{k + k' + 1) \ 
{k,l) £U,ik',l') G V}. 

Theorem 10 f [T7l Theorem 5.4]). Let X be a nonempty finite subset of Ll{d) and let S be a 
lower set. Then the following are equivalent: 

(1) X is an S-code and an S * S-design. 

(2) X is a tight S-code. 

Let X have an angle set D{X) = {ai, ..., a.;}, and set oq = 1- The number s is said to 
be degree. For 0 < i < s, define the binary relation Ri as the set of pairs {x,y) G X x X such 
that x*y = ai. The following theorem shows that for a finite lower set U and a U *Z^-design 
X, X has the structure of an association scheme when the size of is large enough compared 
with s. 


Theorem 11 1 [T71 Theorem 6.1]). Let X be alA *U-design with degree s for a lower setlA. 
Then the following hold. 

(1) \U\ < s + 1. 

(2) If \U\ > s, then X with binary relations defined from angles is an association scheme. 

(3) If \U\ = s 1, then X is a tight U-code. 

To characterize the tight 3-codes, we use the theory of association schemes. We review 
the theory, see more details [2]. 

Let A be a nonempty finite set and let Ri be a nonempty binary relation on A for 
0 < f < s. The adjacency matrix Ai of relation Ri is dehned to be the (0, l)-matrix whose 
rows and columns are indexed by A such that {Afixy = 1 if {x,y) G Ri and {Afixy = 0 
otherwise. A pair (A, is a commutative association scheme, or simply an association 

scheme if the following five conditions hold: 

(1) Aq is the identity matrix. 

(2) Yli=o ~ where J is the all-one matrix. 

(3) = Ai! for some G {0,1,... , s}. 
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(4) AiAj = Yl=oP’i,jAk for i,j G {0, 1,... , s}. 

(5) AiAj = AjAi for any i,j. 

The algebra A generated by all adjacency matrices Ao,Ai ,... over C is called the adja¬ 
cency algebra or Bose-Mesner algebra. 

Since the adjacency algebra is semisimple and commutative, there exists a unique set of 
primitive idempotents of the adjacency algebra, which is denoted by {Eq, i?i,..., Eg}. Since 
{Eq , Ef ,..., EJ} forms also the set of primitive idempotents, we define i by the index such 
that EE = Ef for 0 < i < s. The adjacency algebra is closed under the entry wise product 
o, so we can define structure constants, the Krein parameters q^j, for Eq, Ei,..., Eg under 
entry wise product: 

' ' k=0 

Both sets of matrices {Aq; Ai,..., and {Eq., Ei,..., Eg} are bases for the adjacency alge¬ 
bra. Therefore there exist change of basis matrices P and Q defined as follows; 

S 1 ^ 

Ai = PjiEj, Ej = QijAj. 

j=0 1=0 

We call P and Q the eigenmatrix and second eigenmatrix of the scheme respectively. For each 
0 < i < s, ki := Pio and mi := Qiq are called the i-th valency and multiplicity. 

6.1 Tight complex spherical 3-codes 

Let X be a 3-code in with D{X) = {a, a, /?}, where a is an imaginary number and /? is a 
real number. Note that (p{X) is a real s-code with s = 1 or 2. When d = 1, |X| = |(/>(X)| < 5 
with equality if and only if (p{X) is the regular 5-gon [Uj. In this case, X has the following 
angle set {e^Tri/s ; Q < i < 4}, which implies that X has degree 4. Thus \X\ < 4 holds. 
When d >2 we can easily find real numbers a, b, c such that P{x) = axx -|- b{x -|- x) -)- c is an 
annihilator polynomial of X, namely X is an 5-code, where S = {(0,0), (1,0), (0,1), (1,1)}. 
From Theorem 4.2 in m and dim(Harm(A;,/)) = d^ + 2d., we have the following 

upper bound for 3-codes, also see Table 1 in m- 

Theorem 12. Let X be a 3-code in Ll{d). Then 

\x\< h 

^d"^ P 2d if d> 2. 

However, a tight 3-code is rare, shown in the following theorem. 

Theorem 13. Let X be a 3-code in Ll{d) attaining equality in Theorem 1121 Then one of the 
following holds; 

(1) d = 1 and D{X) = {±-v/—1, —1}, 

(2) d = 2 and D{X) = {±^, -1}. 
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Proof. Let X be a 3-code in r2(l). Since D{X) contains a real number, there exist x, y G X 
such that x*y G M. Since x,y G ^1(1), it must hold that y = —x. Let z,w be the remaining 
elements in X. If z ^ {y^—lx, —^J—lx} held, then X would be a 5-code. Thus {z,w} = 
{y/—lx, —\/—lx} and D{X) is the desired form. 

Let d be an integer at least 2. Since X is a tight 5-code, X is a tight design with respect 
to 5 by Theorem [TOl Since the degree is 3, X carries a commutative association scheme by 
Theoremlllh2i. Then the Gram matrix of X is a scalar multiple of some primitive idempotent 
of the scheme, say Ei. And we arrange the ordering of the primitive idempotents so that 
E 2 = EJ holds and E 3 , is a real matrix. Then 6 = 0, i = 2,2 = 1,3 = 3 hold. 

We will determine the Krein matrix Bl = {qf j)o<j,k <3 the second eigenmatrix Q. 

By 0 Proposition 3.7 (ii), (hi)], we have = gf g = ^ = q^^^ = ^^3 = 0, q^g = 1 

and (? 3,2 = Since Ei = E 2 = and E^ = holds, 

E 2 = o and E^ = o (jXjLii), where Hj^ i is the (/c,/)-characteristic 

matrix of X (see m Section 3] for the (fc, Z)-characteristic matrix) and o (|X|iiii) denotes 
the matrix obtained by applying gk^i to each entry of jXji^i and (see [T71 Section 2 ] for the 
polynomial gk,i{x)). By the recurrence of the Jacobi polynomial (see [T71 Section 2 ]), the 
Krein parameters 912 ) 912)^12 are same as the coefficients of the Jacobi polynomials in the 
product yi^o(^)fl'o,i(3;), namely 9^2 = 2 ~ ^ and 9 i 2 = holds. Since X is a 5 * 5-design 

and 5*5 contains ( 2 ,1), 3 = 0 holds by [171 Corollary 9.3 (ii)]. Since Yl^=o^ij ~ ^ holds 

for any 0 < /c < 3 [21 Proposition 3.7 (v)], we have 

Qi,i + Qi, 3 = d, ( 5 ) 

+ = ( 6 ) 

By the commutativity of the association scheme, q^j = q^^. We have mi = dim(Harm(l, 0)) = 

d and m 3 = dim(Harm(l, 1)) = d^ — 1. Substituting the values mi, m 3 into rrikqfj = 
for {i,j,k) = (1,1,3) [21 Proposition 3.7 (vi)], we have 

(d^ - 1 ) 9?,1 = dql^. (7) 


Using the qf^q{^ 

have 



(9?,i)^ + ^9?,i9p3 


(1,1, 2,1) ([21 Proposition 3.7 (vii)]), 
_ 2 ^ 

d+i ■ 


we 


( 8 ) 


We solve the equations (l5])-([7|) to obtain 


(9i,i) 9i,i) 9 i,3) 91 , 3 ) — 

{ fd{d-{d-l)Vd+2) d^{d+l+VdX2) d{d-l){d+l+Vd+2) d‘^{d‘^-2-Vd+2) ^ 

V dP+dXi ’ {d+l){d:^+d-l)' P+d-l (d+l)(d^+d-l) 

, d{d+{d-l)Pd+2) d?{d+l-^/d+2) d{d-l){d+l-Pd+2) d'^{d'^-2+Pd+2) s 
^ d'^+d-1 ’ (d+T)(d^TTTJ’ P+d-1 ’ {d+l){d'^+d-l) v 


(9) 


First we consider the former case in Q. The Krein number qf ^ is nonnegative only when 
d = 2. When d = 2 holds, the second eigenmatrix Q is given by 


/I 

1 

1 

Vi 


2 2 3 \ 

2 x/^ 2 V^ 1 

V3 V3 ^ 

2 v^ 1 

"Us V3 ^ 
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Next, in the latter case in (l9|), we set t = The second eigenmatrix is given by 


Q = 


/I 

1 

1 


t ^-2 

t+1 


2{t''i-2t+l) 

-6-3t+3t'^+2t^-W-3t'^-2t+5 


t ^-2 

t+1 


-6-3t+3t'^+2t^-tV-3t^-2t+5 

{t^-2){t'^+t-l+W-3fi-2t+5) 
2{t''i-2t+l) 


(f+l)(f2_3) 

t'^+t-l 

(t+l)(f2_3) 

p+t-1 


Then the valency corresponding to the second row of the second eigenmatrix is ki = —— ■ 

By substituting t = ^/d + 2 the valency ki is equal to +6d-5jh4(d-i)y/d+2) ^ which im¬ 

plies that t = y/d + 2 must be an integer. The partial fractional decomposition 243fei = 
81t^ -|- 108t^ — 180^2 — 348t — 149 -|- 3 ^:^ shows that 3t -|- 5 divides 16. Since t is positive, we 
have t = 1 and then d = —1. This contradicts to the fact that d is positive. □ 


For d = 1 , 2 , the tight 3-code is unique, that is proved in Sectional The tight 3-code 
in n(l) is X = {±1, ±-v/—1}. The tight 3-code in 0(2) is {±xi, ±X 2 , ±^ 3 , ±* 4 }, where 
XI = (1, 0), X2 = 1/V6(v^, 1 + V^), X3 = 1/VG{V^, 1 - v^), X4 = 1/VQ{V^, -2). 

Remark 5. For 5 = {(0,0), (1,0), (0,1), (1,1)}, the tight 5-codes with degree 4 were given in 
[m Example 10.2]. They are obtained from the subconstituents of SIC-POVMs in dimension 
d = 2,8. SIC-POVMs are the tight projective 1-codes, see [16] more details. 
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